Non-Markovian quantum trajectories: an exact result 
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We analyze the non-Markovian stochastic Schrodinger equation describing a particle subject to 
spontaneous collapses in space (in the language of collapse models), or subject to a continuous 
measurement of its position (in the language of continuous quantum measurement). For the first 
time, we give the explicit general solution for the free particle case (H = p 2 /2m), and discuss the 
main properties. We analyze the case of an exponential correlation function for the noise, giving a 
quantitative description of the dynamics and of its dependence on the correlation time. 
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The theory of non-Markovian quantum dynamics is a 
subject of growing interest, both from the theoretical 
point of view, as well as for its experimental implica- 
tions On the more theoretical side, interest ranges 
from the theory of open quantum systems to the 

theory of continuous quantum measurement [5[ , to quan- 



to models of spontaneous wave 
With particular reference 



turn trajectories I 
function collapse [ 

to the latter, recent investigations [13| have shown that 
they might need to be generalized to non- white noises, 
in order to be compatible with current experimental con- 
straints. 

Unlike the theory of Markov dynamics in Hilbert 
space, which has been deeply investigated and well un- 
derstood, the theory of non-Markovian quantum dynam- 
ics is still under construction. Important results have 
been already obtained [!, 0] . With particular reference 
to stochastic Schrodinger equations (SSEs) in Hilbert 
spaces, these have been formally generalized to non- 
Markovian noises 0], but explicit results have been ob- 
tained only for simple systems 0, or through approxi- 
mation schemes [l2|, |l4j . 

In the Markovian case, among all SSEs, the following 
equation, 



A 



-Hdt + VXqdW t - -q 2 dt 



(1) 



has received considerable attention If 



17, m 



IS ED, [22 1 . q is the position operator of the particle, H 
its quantum Hamiltonian, Wt a standard Wiener process 
defined on a probability space (QfJ 7 , Q) and A is a pos- 
itive coupling constant [23(. The evolution described by 
Eq. ([1]) is manifestly non-unitary, the non- Schrodinger 
terms being devised in order to reproduce the collapse of 



the wave function [24 1 



The reason why Eq. (JTJ) is so popular is that it rep- 
resents an excellent compromise between mathematical 
simplicity and physical adequacy. From the mathemat- 
ical point of view, it is simple enough to be analyzed 
in great detail [19, 0, HH- From the physical point of 
view instead, it represents a very realistic model describ- 
ing a quantum particle subject to spontaneous collapses 



in space (within collapse models [15|, l2fj), or a particle 



whose position is continuously measured by an appro- 
priate device (within the theory of continuous quantum 
measurement 3, l3|)i or a particle coupled to an envi- 
ronment via its position (within the theory of open quan- 
tum systems [18|). In all fields of applicability, Eq. (Q]) 
has been used to get a deep insight into the dynamics of 
more complicated physical situations. 

It is then of primary interest to study the generaliza- 
tion of Eq. fTJ) to the non-Markovian case. Such a gen- 
eralization has been first proposed in [7j], and reads: 
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VXqwt — 2vAo / dsa(t,s)- — 
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(2) 

where now u>t is a Gaussian non-white noise defined on 
(f2, T, Q), having zero average and the correlation func- 
tion a(t, s). The non-Markovian character is clearly dis- 
played by the third term, which depends on the whole 
past history. For this reason, technically speaking, the 
integration should begin at s = —00. Here we are mak- 
ing the assumption, which actually is an approximation, 
that the state of the system at time suffices to unfold 
the subsequent evolution. This would be the case, e.g., 
if the system has reached an equilibrium configuration 
which is independent of the way it has been reached, and 
one is interested in studying what happens if at time 
the system is driven away from it by a sudden interaction. 

Setting Eq. (T2]) has represented a very important 
achievement. However it remains still somewhat formal, 
as no explicit solutions are known. In this Letter we 
present a recent result, whose technical details are re- 
ported in [2^|: for the first time, the explicit expression of 
the Green's function associated to Eq. © has been com- 
puted, in the case of a free particle (H = p 2 /2m), and its 
properties have been analyzed in detail. The technique 
which has been used can be straightforwardly generalized 
to include linear and quadratic potentials (thus bounded 
systems can also be studied). More complicated situa- 
tions can be analyzed through a perturbation expansion 
on VX 

The Green's function. In [7[ it was first shown that 



2 



the Green's function G(x,t; xq,0) associated to Eq. ([2]) 
allows for the following path-integral representation: 



rq(t)=x 

G(x,t;x ,0) = / V[q]e 

Jq(0)=x o 



S[q] 



(3) 



where the 'action' S[q], which is not standard, having 
both a real and an imaginary part, is: 



S[q] = / ds 



^ q' s 2 + V\q s w s - Xq s J dr a(s, r)q r 



(4) 

We have computed the path-integral in ((3]) using the 
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polygonal approach of Feynman 
is long, in particular due to the last term which contains 
a double integration reflecting the non-Markovian char- 
acter of the evolution; nevertheless the computation can 
be carried out exactly. We report on the final result, fo- 
cusing on the case of a time-translation invariant noise 
{a(t, s) = a(\t — s|)), which is sufficient for most physical 
purposes. In this case, the Green's function becomes [25| : 



■ exp [-A t (xl + x 2 ) + B t x x + C t x + V t x + £ t ] . (5) 

The first two coefficients At and Bt are deterministic 
functions of time and arc defined as follows: 



Ai 
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ti(t), 



(6) 



while the remaining coefficients Ct, T> t and £t depend also 
on the noise Wt through the expressions: 
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St 



^h'M + ^^diwtMi), 



f^(*) + ir/ diwim-i), 



dl wih t (l). 



(7) 
(8) 
(9) 



(Here above and in the following, the symbol ' denotes 
differentiation with respect to the variable within paren- 
thesis.) The random function ht(s) satisfies the following 
non-homogeneous integro-differential equation: 



(10) 



with boundary conditions h t (0) = h t (t) — 0. The func- 
tion ft (s) instead satisfies the homogeneous equation as- 
sociated to Eq. (flO)) . with boundary conditions /t(0) = f , 
ft(t) = 0. Also the function u(t) can be given an an- 
alytic expression in terms of the solution of an integro- 
differential equation; since however the whole square root 
in {5J represents a global factor whose real part looses im- 
portance when normalizing the wave function, and whose 



imaginary part represents an uninteresting global phase 
factor, we omit to write the explicit expression of u(t). 

Eqs. (|5j)- (|I0p represent our main result, from which 
the subsequent discussion follows. One should notice the 
quite remarkable fact that we have been able to compute 
the Green's function associated to Eq. ([2|) (which can be 
applied to any C 2 initial state, giving its time evolution), 
while in general non-Markovian dynamics do not allow 
for such a thing as the Green's function. This fact is less 
surprising if one looks back at how Eq. ([2]) was derived 0] : 
first the evolution was set by means of a propagator, and 
only afterwards the associated differential equation was 
deduced. 

Another relevant observation to make is that the struc- 
ture of the non-Markovian Green's function G(x,t; xn,0) 
is the same as the corresponding Markovian one [2l|, |28( , 
and of course reduces to it in the white-noise limit, as 
proven in [25|. In particular, the exponent is quadratic in 
the variables xq, x, and the coefficients associated with 
the quadratic terms do not depend on the noise. This 
fact has two important consequences: first, the shape of 
Gaussian states is preserved during the evolution; sec- 
ond, their spread evolves deterministically in time. We 
will come back on these points later. Since more gen- 
eral states can be written as superpositions of Gaussian 
states, these facts suggest that any reasonable initial 
state converges almost surely to a Gaussian state with 
a fixed spread both in position as well as in momentum. 
This property holds in the Markovian case, and has been 
subject of an intense investigation [l7, 18, 2(1 21, 2^, 29| . 
It would be important to check it also in a non-Markovian 
setting. 

As a second relevant consequence of Eq. ([5]), one can 
verify that the following ansatz: 
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■4> t = [qa t (s) +pb t (s) +ct(s)]cf>t, 



(11) 



first proposed in [8(, is correct; the three coefficients have 
the following time dependence 251 ]: 



o*(s) = ft(t- 



c t (s) = h t (s) 



Ms) 

2/t'(*) 



Kit) 
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(12) 



(13) 

One can then replace the functional derivative appearing 
in @ with pT|) , giving the non-Markovian equation a less 
cumbersome expression. The form (1111) for the functional 
derivative should make it clear that the non-Markovian 
term of Eq. ([2]) depends on the interplay between the 
Hamiltonian and the collapse terms, since a term pro- 
portional to p appears, which can come only from the 
free part of the evolution. This is the ultimate reason 
why the functional derivative can be computed explicitly 
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only when all o per ators appearing in Eq. Q commute 



a per ai 

■m 



h t (s) — ht(s) — h\ (t) ft(t — s) , where h^(s) is a particular 



or in simple enough cases like ours. solution of (I10|) , namely: 



that the mean position 
l(p)t] evolve accord- 



with each other ^ 

One can further prove [25 1 
JEq [(<?)*] an d mean momentum Ej; 
ing to the classical laws. Moreover, the fluctuations of 
the position of the particle around the average, measured 

by Y q := \Je<q [(q)t — ^[(g)*]] 2 , scale with the inverse 
square root of its mass; this means that, the bigger the 
system, the less random the motion within a given time 
interval. 

Exponential correlation function. The explicit form of 
the coefficients At~£t defining the Green's function de- 
pend on the solution ht(s) of Eq. (fT0|) and on the solu- 
tion ft(s) of the corresponding homogeneous equation. In 
general, this equation cannot be solved explicitly, though 
a perturbation expansion is always possible, which gives 
meaningful results to first orders in A. Nevertheless, the 
solution can be found for particular types of correlation 
functions [3fJ. Among these, the physically most mean- 
ingful example is the exponential correlation function: 



a(t,s) = ( 7 /2)e- 



-~f\t— s\ 



(14) 



where 7 is the inverse of the correlation time. 

With this choice for a(t,s), the homogeneous equa- 
tion for ft(s) can be solved as follows. By differentiat- 
ing twice Eq. (|10|) with w s = 0, one can transform the 
integro-differential equation into the fourth-order differ- 
ential equation [251 ]: 



r{s)-i 2 f"{s)+i^f{s) = i), 



(15) 



where u) = 2^/ftA/m. The general solution is /t(s) = 

Z)fe=i[/t,fc sinh VkS + g t ,k coshu fc s], where f t ,k, gt,k are de- 
termined by the boundary conditions, and v\, v-i are the 
two non-symmetric roots of the bi-quadratic characteris- 
tic polynomial associated to Eq. (fT5|) : 



7,1,2 = V( 7 2 ± 0/2, 



(16) 



Two boundary conditions are already given: /t(0) = 1 
and ft(t ) = 1. The other two conditions can be recov- 
ered [30] from the procedure which led to Eq. (Ti"5"|) and 
read: f»'(0) = 7 / t "(0) and f' t "(t) = - 7 / t "(i). Inserting 
these conditions, one obtains: 



/*(*) 



J2k [ r t sinhu fc (t - s) + U} coshu fc (f — s) — u*] 



^2 k [2c -I- r\ sinhufct + u\ coshufct] 



(17) 

with k = 1,2 and where r\ = aj.cosh.v^t + bj.smhv^.t 
and Uf = dkSmhvfct — ccoshv^t; we have also defined: 
ak = 7^[^ + (-l) fe C], b k = ^[«J+(-l)Va c = 
4 = -JVk v h with = 2 if A; = 1, fc = 1 if fc = 2. 

The function ht(s) can be found in a similar way, 
though its expression is more complicated, as h t (s) solves 
the whole inhomogeneous equation. Taking into ac- 
count the boundary conditions, h t (s) takes the form: 



Ml) 



Ul) {w'{ - ^ Wl ) dL 



sinhui(s — I) svah-Vzis — I) 



(18) 



The problem has been completely solved. One can check 
that in the white-noise limit 7 — > 00 (a(t, s) — > 5(t — s)), 
one recovers the well-known Markovian expressions. 

Evolution of Gaussian states. The analysis of Gaus- 
sian states is particularly useful in order to understand 
the behavior of a typical physical state. As previously 
anticipated, the shape of Gaussian wave functions does 
not change in time. In fact, an initial state: 



4>o{x) = exp[-a a; 2 + fox + 70] 



(19) 



preserves its functional dependence on x, while the com- 
plex parameters cto, fo and 7 o evolve in time as follows: 



a 1 



It 



At 



-fo + £ t + 



(C t + A>) 2 



Ct+Po 
'4(ao +-4*) 



A(a + At 



(20) 



Analyzing the above expressions with the help of 
Eqs. ©-([5]), one immediately sees that the evolution of 
at is deterministic, while (3 t and 74 have stochastic terms. 
This means that, like in the white-noise case, both the 
spread in position and in momentum of <j) t (x) , which are 
given by a t , evolve deterministically in time. On the 
other hand, both the mean position and the mean mo- 
mentum, which depend both on a t and fit , have stochas- 
tic components; their stochastic averages instead evolve 
according to classical laws, as we have already antici- 
pated. 

We focus now our attention on the spread in position 



a(t) = 1/2, 



in the case of the exponential correla- 



tion function treated before. Fig.[T]shows how the spread 
evolves, for different values of 7. Qualitatively the behav- 
ior is the same for any 7: the wave function shrinks in 
space, reaching an asymptotic finite value. On a more 
quantitative level, we see that the stronger 7, the faster 
the collapse. One can also notice that the collapse is ef- 
fective starting with relatively small values of 7: a value 
7 ~ 10 sec -1 already ensures that after about 10~ 3 sec 
the wave function has collapsed below 10™ 5 cm, which is 
the threshold chosen by GRW 24J, below which a state 
can be considered as localized. This means that the pos- 
sibility opens for non-Markovian models to be as effec- 
tive as the corresponding white-noise models as far as 
the collapse process is concerned, but, at the same time, 
to give different physical predictions regarding specific 
experimental situations. This possibility has first been 
suggested in [l3|. 
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FIG. 1: Time evolution, under the assumption of an expo- 
nential correlation function, of the spread in position a(t) of 
a Gaussian wave function. <r(0) has been set = 1 m. The 
value 7 = oo corresponds to the Markovian case. The other 
parameters have been chosen as follows: m = 1 Kg, Ao = 10~ 2 
m -2 sec -1 . Time is measured is sec, distances in m. 



From the previous expressions one can explicitly com- 
pute the asymptotic value of at, which is: 

Z7TL 

cico = lim a t = - — {vx + v 2 -j). (21) 
t^oo zn 

The quantity l/2^/a^ is the final spread in position to 
which all Gaussian states (and, reasonably, any initial 
state) converge to, in the long-time limit. 

Conclusion. We have computed for the first time the 
Green's function associated to the motion of a free par- 
ticle as described by Eq. ([2]), from which the entire non- 
Markovian dynamics can be unfolded. We have analyzed 
the physically important case of an exponential correla- 
tion function. By studying Gaussian states, we have seen 
how the collapse occurs, and have derived an exact ex- 
pression for the asymptotic spread. The tools we have 
employed to derive the above results are flexible and can 
be applied to more complex physical situations. 
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